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ABSTRACT 

A generic feature of weakly interacting massive particle (WIMP) dark mat- 
ter models is the emission of photons over a broad energy band resulting from 
the stable yields of dark matter pair annihilation. Inverse Compton scattering 
off cosmic microwave background photons of energetic electrons and positrons 
produced in dark matter annihilation is expected to produce significant diffuse 
X-ray emission. Dwarf galaxies are ideal targets for this type of dark matter 
search technique, being nearby, dark matter dominated systems free of any astro- 
physical diffuse X-ray background. In this paper, we present the first systematic 
study of X-ray observations of local dwarf galaxies aimed at the search for WIMP 
dark matter. We outline the optimal energy and angular ranges for current tele- 
scopes, and analyze the systematic uncertainties connected to electron/positron 
diffusion. We do not observe any significant X-ray excess, and translate this 
null result into limits on the mass and pair annihilation cross section for parti- 
cle dark matter. Our results indicate that X-ray observations of dwarf galaxies 
currently constrain dark matter models at the same level or even more strongly 
than gamma-ray observations of the same systems, although at the expenses of 
introducing additional assumptions and related uncertainties in the modeling of 
diffusion and energy loss processes. The limits we find constrain portions of the 
super symmetric parameter space, particularly if the effect of dark matter sub- 
structures is included. Finally, we comment on the role of future X-ray satellites 
(e.g. Constellation-X, XEUS) and on their complementarity with GLAST and 
other gamma-ray telescopes in the quest for particle dark matter. 
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Introduction 



The fundamental nature of dark matter is at present unknown. It is widely believed that 
dark matter is in the form of a particle lying outside the ranks of the Standard Model of par- 
ticle physics. An attractive possibility is that the New Physics sector hosting the dark matter 
particle is connected to the electro-weak scale, soon to be explored with the Large Hadron 
Collider. Motivations in support of this possibility include the fact that several models for 
new, electro-weak scale physics encompass particles that have all the micr oscopic features 



neces sary to be the dark matter (this is the case for su persvmmetrv fsee e.g. lJungman et al. 



1996 ). models w ith universal extra-d imensions (see e.g. Hooper fc Profumo 2007h . lit tle Higgi 



models (see e.g. iBirkedal et al.ll2006l ). and many others (see e.g. iBertone et al.ll2005l )): in ad- 
dition, weakly interacting particles with electro-weak scale mass (or WIMPs, for weakly 
interacting massive particles) possess the virtue that an order of magnitude estimate of their 
thermal relic abundance falls in the same range as the density of cold dark matter inferred 
from various cosmological observations f Komatsu et al.l boosi ) . WIMPs might have once 
been in thermal equilibrium in the early universe, eventually freezing-out when the rate of 
pair-annihilation fell below the Hubble expansion rate: particle dark matter, in this scenario, 
would be one more ^""thermal relic" , similar to the cosmic microwave background photons or 
the light elements produced during Big Bang Nucleosynthesis. 

The fact that, occasionally, WIMPs can still pair-annihilate in the present universe fos- 
tered a wide array of indirect particle dark matter searches: the pair annihilation of dark 
matter is expected to yield, for instance, a significant amo unt of energetic antimatter (includ- 



ing GeV positrons, antiprotons and antideuterqns (se e e.g. iBaltz fc Edsiolll999l:lDonato et al. 



2004; Profumo fc UUio 2004; Baer fc Profumo 2005; Brun et al 



2007 



Bringmann &: Salati 



20071 ). and gamma rays, from various radiative processes and from the decay of neutral 
pions resulting from the hadronization of particles in the pair-annihilation final state. Of 
interest to us here is the fact that particle dark matter annihilation also produces a popula- 
tion of energetic electrons and positrons from, for instance, charged pion, muon, gauge and 
Higgs boson decays, and, possibly, prompt production. The injection spectrum depends on 
the details of the particle dark matter model, but since electrons and positrons (e"'"e~) are 
produced at a center-of-mass energy corresponding to twice the particle dark matter mass, 
assumed here to be around the electro-weak scale, many of them will have energies at or 
above a GeV. These energetic e~^e~ populate dark matter halos in any generic WIMP model. 
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with densities which depend on both the dark matter density profile and the WIMP pair 
annihilation rate. Electrons (and positrons) diffuse, loose energy and produce secondary 
radiation through various mechanisms. In the presence of magnetic fields they emit at radio 
wavelengths via synchrotron radiation. Inverse Compton (IC) scattering off target cosmic 
microwave background photons and background light at other frequencies gives rise to a 
broad spectrum of photons, stretching from the extreme ultra-violet up to the gamma-ray 
band. A further, typically subdominant contribution to secondary photon emission results 
from non-thermal bremsstrahlung, i.e. the emission of gamma-ray photons in the deflection 
of the charged particles by the electrostatic potential of intervening gas. 

The mult i- w avelength em i ssion f rom dark matter annihilation was st udied in detail in the 



semin al works of iBaltz &: Wail (120041 ) for galactic dark matter clumps, i n IColafrancesco et al. 



( I2OO6I ) for the case of the Coma cluster and in IColafrancesco et al.l (l2007al ) for the dwarf 
spheroidal galaxy Draco. Other recent studies include an interpretation of the significant 
non-thermal X-ray activity o bserved in the O phiuchus cluster in terms of IC scattering of 
dark-matter produced e~^e~ (jProfumol l2008al ). an analysis of the broad-band dark-rn atter 
annihilation spectrum expecte d from the Bullet c luster (IColafrancesco et al.l l2007bl ) and 
from the galactic center region (IRegis fc Ullioll2008l ). In addition, radio emission from e+e~ 
produced in dark matter annih ilation was c onsid ered as a possible source for the 'WMAP 
h aze" in the semin al paper of iFiiikbeiner J2004h. and subse quently analyzed in detail in 



Hooper et al.l (120071 ) , iHooper et al.l (120081 ) and iHoopen (120081 ) . Other studies have also pre 



viously addressed synchrotron radiation indu ced by dark matter annihilation (e.g. 
2OO0I : iBertone et ZII2OO1I : lAloisio et al]l2004l ). 



Gondolo 



Among the possible targets for the observation of an astronomical signature of dark mat- 
ter annihilation, local dwarf spheroidal (dSph) galaxies stand out as excellent candidates for 
several reasons. First, unlike the galactic center region or galaxy clusters, no significant 
diffuse radio. X-ray or gamma-ray emission is expected: the gravitational potential well of 
dSph galaxies is too shallow for them to host any sizable thermal bremsstrahlung emission 
at X-ray frequencies, and, more importantly, the gas densities appear to be extremely low 



(see e.g. lMatedll998l ). Second, dSph are the most dark matter dominated known systems, 
and, with the exception of our own Milky Way, they are the closest known bound dark 
matter systems. Unlike a signal from the galactic center region or from a nearby cluster, a 
diffuse X-ray or radio emission from a nearby dSph galaxy would likely not have an obvious 
astrophysical counterpart that could fake a dark-matter induced emission. The cross corre- 
lation of diffuse emission from one of the nearby dSph galaxies with, for instance, point-like 
emission at gamma-ray frequencies detected with GLAST (unlike secondary emission from 
e~^e~, species which undergo spatial diffusion, gamma rays trace the dark matter density 
profile squared), and a study of the spectral features at various wavelength could potentially 
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lead not only to the conclusive detection of dark matter, but even to the identification of 
the particle mass and some of its microscopic particle properties. The spectral shape of the 
broad-band emission depends in fact on the dominant Standard Model final state into which 
the dark matter particle annihilates, and this is in turn determined by the WIMP particle 
model. 

In this paper, we present the first systematic study of archival X-ray data on local dSph 
galaxies with the aim of detecting a signal from dark matter annihilation. After introducing 
the physics of e~^e~ production from dark matter annihilation and of subsequent diffusion and 
energy loss, we present in Sec. [2] a few examples of multi-wavelength spectra, emphasizing 
and analyzing the role of systematic uncertainties connected with the dark matter particle 
model and with the diffusion setup. We then present our data reduction and analysis in 
Sec. El motivating our choice for the energy range and angular region. We present limits 
on the X-ray diffuse emission from the dSph galaxies Fornax, Carina and Ursa Minor in 
Sec. 14. 1[ Sec. 14.21 is devoted to a study of how our null result limits particle dark matter 
models. Our summary and a discussion of the role of future X-ray detectors in particle dark 
matter searches, including their connection with gamma-ray telescopes, are presented in the 
final Sec. O 



2. The Multi-wavelength Spectrum from Dark Matter Annihilation 

The starting point for the computation of the multi- wavelength spectrum resulting from 
dark matter annihilation is the identification of the relevant source function. Given a stable 
sp ecies i and a position r, th e source function Qi{r, E) (we are using here the same notation 



of IColafrancesco et al.ll2006[ ) gives the differential number of particles i per unit time, energy 



and volume element produced at r and with an energy E\ 

AAff 

Q^if, E) = {av)o J2 -T^iE) Bf (1) 
/ 

where {(Jv)q is the WIMP pair annihilation rate at zero temperature, and the sum is over 
all kinematically allowed Standard Model annihilation final states / (for instance, quark- 
antiquark, W^W~, lepton-antilepton etc.), each weighed with a branching ratio Bj and 
producing a spectral distribution dNf /dE, after prompt production or decay and fragmen- 
tation into the stable particle i. Finally, A/'pairs(^) is the number of dark matter particle pairs 
per volume element squared, which for the case of a smooth dark matter distribution Pdm(0 
is given by 

A/-p.„(r-) = (2) 
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where we indicate with moM the mass of the dark matter particle. 



In the case of gamma rays, since photons propagate on straight hnes, the flux from 
prompt emission in a given direction is simply given by the integral of the appropriate 
source function along the line of sight, 



dl Q^{E^,f{l)). 



(3) 



l.o.s. 



The quantity in Eq. is then integrated over the angular region over which the signal is 
observed, and convolved with the angular dependent instrumental sensitivity of the gamma- 
ray telescope under consideration. 

The treatment for electrons and positrons is complicated by diffusion and energy loss 
processes. We model these with a diffusion equation of the form 



d due 
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(4) 



where drie/dE is the electron and positron spectrum, D{E,r) is the diffusion coefficient 
and h{E, r) is the energy loss term. Ana lytical solutions to the equ ation above exist in the 
equilibrium regime (see e.g. App. A in IColafrancesco et al.l 120061 ) . as long as the spatial 
dependence of both the diffusion coefficient and the energy loss term are dropped, and 
spherical symmetry in the dark matter distribution is assumed (in the present study we make 
the same hypotheses). The dependence of the diffusion coefficient on energy is assumed to 
be a power law of the form 

Very little is known about diffusion in systems of the type of interest here, i.e. dSph galaxies. 
This forces us to make educated guesses about Dq and 7. The best kn own system as far as 
cosmic ray propagation is concerned is by all means our own Galaxy. iDonato et al.l (120041 ) 
analyzed data on cosmic ray fluxes in the Milky Way, in particular ratios of primary to 
secondary species. The outcome of their analysis was to determine the preferred values for 
the parameters Dq and 7, in the framework of a diffusion setup similar to the one outlined 
here. Their median values are Dq = 1.1 ^ 10^^ cm^/s and 7 = 0.7, which we emplo y here as 
reference values. Using the phenomenological ranges found by iDonato et al.l (120041 ). we will 
also consider values of the parameter < 7 < 1. 

The parameter Dq is related to the size and scale of galactic magnetic field inhomo- 
geneities; data on galaxy clusters indicate that larger systems feature larger values for Dq 
than galactic-size systems. This can be understood by estimating the diffusion coefficient 
Dq in the context of passive advective transport in a turbolent flow as the product of the 
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turbolent velocity Vr, and o f the turbolent injection scale L (see e.g. iNarayan &: Medvedev 
2OOII : iLazarian fc Clidl2004l ). If the turbolent velocity in clusters of galaxies is of the order 
of the velocity dispersion of galaxies (say Vl ~ 1000 km/s), and the injection scale is of the 
order of the galactic scale (say, L ~ 20 kpc), then the predicted Dp ~ LV[^ 6 x 1 0^° cm^/s 



is indeed close in value to what is observed in clusters (IZakamska fc Narayanll2003l . quote for 
instance Dq ~ 6.22 x 10^" cm^/s for the Hydra A galaxy cluster). Scaling this estimate from 
cluster down to g alactic scales, and usi ng the reference value for Dq from galactic cosmic ray 
data obtained by iDonato et al.l (120041 ) indicates, assuming Vl for the galaxy is smaller by 
a factor 10 than in clusters, that the value of L for galaxies is at least a factor 100 smaller 
than for clusters. Conservatively assuming that L does not change switching from a Milky 
Way size galaxy to a dSph, the simple scaling in Vl of the diffusion coefficient points to 
Dq ~ 10^^ cm^/s, given that the velocity dispersion of the Milky Way is more than one order 
of magnitude larger than those observed in local dSph galaxies. 

In this respect, we might expect values for Dq smaller than the conservative Milky Way 
median value we use as a reference here. We will therefore also consider as an alternative 
value Dq = 10^^ cm^/s, and we will discuss the dependence of the X-ray emission from dark 
matter pair-annihilation on both Dq and 7 in Sec J4.2[ 

Another important parameter for the diffusion model is the diffusion volume, where 
Eq.(jl]) is solved, at the boundary of which free-escape boundary conditions are imposed. 
Following IColafrancesco et al.l (l2007al ). we assume a spherical diffusion zone, and a diffusion 
radius Vh corresponding to twice the radius of the stellar component, typically a few kpc 
for local dSph galaxies (jMatedll998l ). This choice is again motivated by analogy with the 
picture for the Milky Way. The value of Vh is also a crucial parameter for the computation 
of the X-ray emission from dark matter annihilation: the smaller the radius of the diffusion 
region, the larger the fraction of e~^e~ escaping and diffusing away, hence the smaller the 
X-ray emission. 

The e nergy loss term, which is a ssumed to be position- independent, includes the follow- 
ing terms (IColafrancesco et al.ll2006l ): 



b{E) = hc{E) + bsyn{E) + 6coul(^) + b 

brem 



(6) 



where 
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n {1 + log(7e/n)/75) 



6?c ^ 0.25 X lO"^*^ GeVs"^ 



6°y„ ~ 0.0254 X 10"^^ GeVs"^ 



~ 6.13 X 10-^^ GeVs"^ 
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6brem(^)=Cm(M7e)+0.36) 
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1.51 X 10"^^ GeVs 



where n indicates the thermal electron density and 7e = E/lmeC^). We set in what follows the 



average magnetic fie 
gala xies reported in 



see 



d toB 



Klein et al. 



1 fi G, in concordance with radio observations at 5 GHz of dSph 

6 



Colafrancesco et al.ll2006l ) 



(I1992I ). and the thermal electron density to n = 10 cm 



No measurements are available on the average magnetic field for the dSph galaxies under 
investigation here. Eq. ([6]) indicates that with our reference choice for B, energy losses for 
energetic electrons and positrons (£" > 1 GeV) are dominated by Inverse Compton processes. 
Smaller values for the magnetic field would thus in no way affect our results. However , 
larger values for B could be in principle allowed by available data (see e.g. iKlein et al.lll992l ). 
Synchrotron losses dominate over IC for i? > 4 /iG. Above that value, we expect a suppression 
on the secondary X-ray emission from WIMP annihilation, since the density of energetic 
electrons and positrons are depleted by more efficient energy losses. Larger average magnetic 
fields would however also produce an intense radio emission, also f rom se condary processes 
related to WIMP annihilation. The results of I Colafrancesco et al.l (l2007al ) indicate that, in 
the case of the Draco dSph, for 5 3> 1 /iG, the synchrotron emission originati ng from WIMP 
annihilation would violate VLA limits on radio emission (IKlein et al.lll992l ). The thermal 
electron density only affects the strength of Coulomb losses (see Eq. (|6])), and thus the low 
energy part (£" -C 1 GeV) of the electron-positron equilibrium spectrum. Larger values for 
n would thus only affect the low energy end of the IC emission in the mult i- wavelength 
spectrum, well below the X-ray regime of interest here. In addition, n ^ 10^^ cm^"^ would 
lead to some significant thermal bremsstrahlung emission, potentially in confiict with the 
limits on X-ray emission we present below. 

After specifying the dark matter density profile and the e^e~ injection spectrum, Eq.Q 
can be integrated to find the equilibrium distribution dne/dE. In turn, knowledge of the 
spatial and energy distribution of e+e~ allows us to compute the multi-wavelength secondary 
emission. For our purposes, and in the energy range of interest here, the dominant contri- 
bution comes from the up-scattering of CMB photons. While a contribution from starlight 
and background light at other frequencies is also expected, it is generically subdominant in 
the X-ray band. 

The inverse Compton power is obtained by folding the number density of target photons 
n{e) (here, the CMB black body spectrum) with the IC scattering cross section: 



P{E^, E) = cE^ J de n{e)a{E^, e, E), 



(7) 



where E = '-yf.rrif.c? is the e+e energy, e is the target photon energy, E^ is the energy 
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of the up-scattered photon, and a indicates the Klein-Nishina formula (the full relativistic 
QED form of the Thomson scattering cross-section). Folding the IC power with the e+e~ 
equilibrium distribution, we get the local emissivity 

,iE-„r)=jAE[^ + '^)p(E,,E). (8) 

Finally, the integrated flux density spectrum, as in the case of gamma rays, is given by the 
angular and line of sight integral of the expression above (see Eq. ([3])). 

The last ingredient needed to actually compute the broad band spectrum of dark matter 
annihilation is to specify the particle dark matter model. Up to a normalization factor 
depending on the dark matter density profile and on the dark matter pair annihilation rate, 
all that matters as far as the particle dark matter model is the mass ttidm and the set 
{S/} of branching ratios into given Standard Model final states. For simplicity and to allow 
comparison with other studies, we choose particle dark matter models with branching ratio 
1 into a certain final state (i.e. always annihilating into the same Standard Model final 
state). We choose as a reference model a 100 GeV WIMP pair-annihilating into a hh pair 
(the resulting spectrum for other quark-antiquark final states is very similar to this one). 
Our choice is also motivated by the fact that many supersymmetric models feature this 
particular final state as the dominant one (examples include the bulk and fu nnel regions of 



minim al supergravity, and generic bino-like models with large tan /3, see e.g. iBertone et al. 



20051). 



We show the spectral energy distribution (SED) for this particular dark matter particle 
model with a black line in both panels in Fig{Tl where we only show photon energies larger 
than 10 eV. We chose the normalization so that the integrated gamma-ray flux between 0.1 
and 10 GeV is 10~* photons per cm^ per s, over an an gular region of one de gree, which roughly 



corresponds to the EGRET point-source sensitivity (IHartman et al.lll9991). For refere nce, we 



use a Navarro- Frenk- White (NEW) dark matter density profile (INavarro et al.l 119971 ) 



\r\\ ^ f\r\ ^ ^ 



PDM(|r1) = p. (^^J V^^7 ' 

where \r\ indicates the distance from the center of the dSph galaxy, and where we set the 
scale radius = 1 kpc. Also, we set Vh = 2.4 kpc. In the left panel we assess how the particle 
physics model affects the dark matter annihilation SED. The dashed green line indicates the 
result for a model with the same 100 GeV mass, but annihilating into t^t~ , a final state 
also motivated by supersymmetry (e.g. in the stau coannihilation region), which features a 
harder e^e~ injection spectrum, as well as a harder gamma-ray spectrum. A more abundant 
population of energetic electrons and positrons results in an IC spectrum which peaks at 
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Photon Energy [GeV] Photon Energy [GeV] 



Fig. 1. — Left: The multi- wavelength spectrum (spectral energy distribution) of various 
particle dark matter models (the reference 100 GeV WIMP annihilating into bb, two setups 
also with a mass of 100 GeV but pair-annihilation final states r+r" and W^W~, and a model 
with the same final state as the reference model, but a mass of 1 TeV). Right: a comparison 
of the effect on the dark matter multi-wavelength spectrum of the reference model of changes 
to parameters in the diffusion model (a diffusion coefficient set to 10^^ cm^/s instead of the 
reference value 1.1 x 10^'' cm^/s, a scaling with energy of the diffusion coefficient 7 = 
instead of the reference value 7 = 0.7, and a radius for the diffusion region of = 0.5 kpc, 
versus the reference value of 2.4 kpc). 
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larger energies than for a softer injection spectrum, such as the reference model. In the soft 
X-ray band, however, while we find different spectral indices, the two models have a similar 
level of emission at a given gamma-ray luminosity. The dot-dashed red line indicates, again 
for = 100 GeV, the case of the W~^W~ final state, found e.g. in Higgsino and Wino 
super symmetric dark matter models. This case is intermediate between the two previous 
ones, and, compared to the reference model has a population of energetic electrons and 
positrons resulting from e.g. W ^ lui, I = e, /i, r, which is responsible for the bump in 
the IC emission in the MeV region. Changing the mass of the dark matter particle shifts 
the SED: a larger yields more energetic e"'"e~, as can be appreciated comparing the 

black and the double-dotted-dashed blue line (moM = 1000 GeV). Again, the X-ray region 
is however only mildly affected. Notice that a generic feature of the X-ray spectrum is a 
hard (0.8 < a < 1.5) photon spectral index, harder than most astrophysical X-ray sources, 
particularly diffuse thermal emission. 

In the right panel of Fig. [1] we study the impact of changing the parameters in the 
diffusion model on the dark matter annihilation SED. A smaller diffusion coefficient (orange 
dashed line, Dq = 10^^ cm^/s) suppresses e+e^ diffusion and escape from the diffusive 
region, yielding a larger IC flux, particularly in the X-ray band of interest here. The choice 
of Dq is therefore very significant in the computation of the expected X-ray flux; we will 
quantitatively analyze this statement in Sec J4.2l FigjSl A milder rigidity dependence in the 
diffusion coefficient (i.e. a smaller 7, set to zero for the double-dashed-dotted indigo line) 
leads, at fixed Dq, to a suppressed diffusion of the most energetic e~^e~ (i.e. the diffusion 
coefficient is smaller for > 1 GeV). As a consequence, we get an enhancement of the high- 
energy IC photons and a suppression at lower up-scattered photon energies. The result on 
the X-ray emission is a mild suppression. Finally, the magenta dot-dashed line indicates the 
effect of taking a diffusive region with a radius a factor 5 smaller than our benchmark choice. 
This corresponds to a diffusion volume more than a factor 100 smaller: in turn, this implies 
a much larger loss of e~^e~ escaping the diffusion region, leading to a significant suppression 
of the electron/positron number density and, therefore, of the IC X-ray signal. 

In summary, significant uncertainties affect the computation of the IC X-ray emission 
resulting from dark matter annihilation; while for a given diffusion model the differences in 
the soft X-ray band are rather mild, changing the parameters in the diffusion setup affects 
quite significantly both the normalization and (although less dramatically) the spectrum 
of the predicted SED; the reference diffusion setup chosen here gives rather conservative 
predictions for the X-ray flux; the outstanding features of the signal we are after are therefore 
(1) extended emission and (2) a hard spectral index. 
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3. Data and Data Reduction 

Of the nearby Local Group dwarf spheroidal (dSph) galaxies, three have observations 
in the XMM-Newton archive: Ursa Minor observed for 52 ks between August and October 
2005 (ObsIDs 0301690201, 0301690301, 0301690401, 0301690501), Fornax observed for 104 
ks in August 2005 (ObsID 0302500101), and Carina observed for 42 ks in May 2004 (ObsID 
0200500201). Here we focus on XMM observations as its large effective area and large field- 
of-view are ideal for detecting faint diffuse emission, while its fairly small PSF allows us 
to exclude X-ray point source contamination from X-ray binaries and background AGN. A 
couple of additional local dSphs have been observed with Chandra and Suzaku, but much 
of these data are not yet public. All of the XMM observations were taken in Full Frame 
mode; for Ursa Minor and Carina the thin optical blocking filter was used, while for Fornax 
the medium filter was used. Unfortunately, the background flare filtering, discussed below, 
reveals that three of the observations of Ursa Minor are highly contaminated with background 
flares. In our analysis, we use only ObsID 0301690401 for this dwarf. In addition, for the 
Fornax and Ursa Minor observations CCD6 on MOSl was not available, but this CCD does 
not fall within the considered source region (see below). 

The data were reduced using XMMSAS version 7.1.0, and all observations were repro- 
cessed using the EPIC chain tasks. For EMOS data, we use patterns 0-12 and apply the 
^XMMEA_EM flag filtering, and for EPN data, we use pattern s 0-4 and flag equal to zero. 



Due to the time variability in the spectra of background flares (iNevalainen et al.ll2005l ). we 
filter for periods of high background in several energy bands. We first apply a conservative 
cut on the high energy count rate matching the 2XMMp pipeline and the EPIC blank-sky 
event file filtering (< 2 cts s~^ for EMOS data and < 60 cts s~^ for EPN data). This cut 
removes the most egregious flares. In the case of Carina, we found that the observations 
contained a number of milder flares that were not removed by this conservative cut, but 
never the less biased the background rate high, and for this dwarf we found it necessary to 
apply more stringent high energy count rate cuts of < 0.65 cts s~^ for EMOS data and < 3.5 
cts s~^ for EPN data. For all dwarfs, we then applied a 3cr clipping to the source-free count 
rate in three energy bands, 0.5-2 keV, 2-5 keV, and 5-8 keV. Here time bins (bin size of 100 
sees) with rates more than 3cr away from the mean are removed recursively until the mean 
is stable. As noted above, the flare filtering excluded nearly all of the exposure time from 
three of the four Ursa Minor observations. The final clean exposure times are listed in Table 
1 along with the adopted dwarf central positions. 
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4. Results 

We describe below the limits we obtain for the diffuse X-ray emission from the se- 
lected dSph galaxies (Sec. 14.11) . and how these limits constrain particle dark matter models 
(SecIO}. 



4.1. X-Ray Flux Limits 

Before comparing our model predictions to the actual X-ray data, we preformed a de- 
tailed study of the optimal energy and angular range where we expect the largest possible 
signal-to-noise ratio (S/N). Specifically, considering both the XMM- Newton affective area 
versus energj0 and the range of possible X-ray spectra resulting from the secondary emissions 
of e^e~ produced in WIMP annihilations, we find that the expected signal-to-noise peaks in 
roughly the 0.5-8 keV band. We illustrate this point in Fig. [2], where we show contours of 
constant signal-to-noise in the plane defined by the lower energy limit Ei and by the upper 
energy limit E2. We integrate over the [Ei,E'2\ interval the signal obtained from WIMP 
annihilations (after diffusion, modeled according to the reference setup described above) and 
we divide by the square root of the background over the sa me energy interval. W e show the 



XMM background count rate in the inset of the left panel (ICarter fc Read! 120071 ). The S/N 
grows towards the red cross, which indicates the best values. In the left panel we assume 
a soft e^e~ injection spectrum (for reference, a 100 GeV WIMP pair annihilating into hh 
pairs), while in the right panel we take a hard injection spectrum (100 GeV WIMP annihi- 
lating into r+r^ pairs). We verified that varying the WIMP mass doesn't affect the position 
of the best energy interval. Also, any combination of final states in the context of particle 
dark matter models motivated by beyond the Standard Model physics falls in between the 
two final states under consideration here. We obtain that with a soft injection spectrum the 
best energy range (giving the highest signal-over background) is [0.44,7.8] (all energies are 
in keV), while for a hard spectrum it is [0.83,9.0]. We decided to use the 0.5-8 keV range as 
the optimal X-ray band for XMM. 

Similarly, we use the XMM sens itivity versus off-axis angle and th e observed dark matter 



profiles of the dwarf galaxies (see e.g. lStrigari et al.ll2007al : lMatedll998l ) to explore the optimal 



source radius. The total size of the dwarfs extends beyond the XMM field-of-view, but the 
expected X-ray fiux decreases significantly with radius as the dark matter density drops as 
well as through the telescope vignetting. Neglecting the effect of diffusion, the maximal S/N 



http://xmm.esac.esa.int/external/xmm_user _support/documentation/uhb_2.5/index. html 
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Fig. 2. — Curves at constant signal over the square root of the background in the plane 
defined by the upper and lower limits Ei and E2 of the energy interval over which the signal 
is integrated. The inset in the left panel shows the average XMM background count level 
we used. The left panel assumes the reference model (100 GeV, bb final state), while the 
right panel uses a model with a harder e~^e~ injection spectrum. The red crosses indicate 
the [Ei,E2] interval giving the maximal S/N. 
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is obtained for the smallest possible angular regions for profiles which diverge with — > 0, 
such as the NFW profile. For cored profiles the best S/N (again in the limit of no spatial 
diffusion) is achieved around the profile scaling radius r^. Factoring in diffusion, we find that 
the optimal angular region significantly increases. The optimal radius therefore depends 
both on the assumed dark matter profile and on the diffusion model, but in all cases we find 
that a radius of ~ 6' is a good choice. 

We investigate whether we detect diffuse X-ray emission from these dwarfs above what 
we expect from the X-ray background. First, we detect and exclude X-ray point sources; 
point sources are detected using the SAS wavelet detection routine, ewavelet, in a mosaiced 
0.5-8 keV band image of all three detectors, and data within a 25" radius region of the source 
locations are excluded. We then compare the detected flux within a 6' radius aperture of 
the dwarf center to the expected X-ray back ground flux in that r egion of the detector by 



comparing to the EPIC blank-sky event files (jCarter fc ReadI 120071 ). The blank-sky files are 



collected from all over the sky, while the dwarf galaxies considered here all lie outside of the 
galactic plane in regions of fairly low hydrogen column density. We, therefore, use the tool 
BGSelector to create blank-sky files using only regions with similar galactic hydrogen column 
density to the dwarf galaxies; specifically, we filter on nH between 10^° and 5 x 10^° cm~^. 
We then apply the same pattern, fiag, and multi-band fiare filtering to the background 
files as was used for the dwarf galaxy observations (Sec. 3). Finally, we re-project the 
appropriate blank-field event files (medium filter for Fornax and thin filter for Ursa Minor 
and Carina) to match the sky position of the dwarf galaxy observations using the routine 
skycast. Renormalizing the blank-field count rate in the source region (i.e. within r = 6' 
excluding point source regions) using the ratio of the count rate in the blank-fields to that 
in the dwarf galaxy observations in an outer region of the detector (where the source fiux is 
expected to be much lower), we find that none of the three dwarfs show significant diffuse 
X-ray emission above what is expected from the X-ray background. 

As shown below, the non-detection of diffuse X-ray emission from dwarf galaxies places 
limits on the possible particle dark matter model. We place an upper limit on the possible 
diffuse X-ray emission from each dwarf galaxy by assuming that all of the detected flux in 
the source region stems from the X-ray background and determining the necessary flux for a 
diffuse source to be detected at 3a above this background. The derived flux limits are listed 
in Table 1 and range between 1 x 10~^ and 2.5 x 10^^ photons cm~^ s~^. 
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4.2. Constraints on Dark Matter Models 



We summarize and compare the diffuse X-ray ffux limits we obtain for Fornax, Carina 
and Ursa Minor in the left panel of FigjSl To model the dark matter density d istribution for 



these t hree dSph galaxies we employ NEW profiles, and follow the analysis of IStrigari et al. 



( l2007al ) for the ranges of scale radii and densities allowed by dynamical data and CDM 
structure formation theoretical constraints. 

Tab. [2] collects the reference, minimal and maximal values for the scaling density p^, the 
scaling radius r^, as well as the reference distance. Tab. [3] indicates instead the reference, 
minimal and maximal values for the angle-averaged line-of-sight integral J for a solid angle 
Ar2 10~^ sr, corresponding to an angle 9 = 6'. The quantity J is defined as 



J 



An 



pDM[r{s)] ds. 



(10) 



l.o.s. 



We express J in units of 10 GeV cm~ . Although we do consider the variation of the 
detector sensitivity with the offset angle (see Sec. [3]), the values of J give an idea of the 
normalization of the signal from different dwarfs (see also IStrigari et al.l l2007al . for further 
details). For instance, given the number A^^ of photons produced per DM annihilation in 
a given energy band and for a certain particle dark matter model, the expected flux 0^ in 
photons cm~^ s~^ in that energy band and from the solid angle AQ from a dSph with a 
normalization J will simply be given by 



0^ = K 



^ "'dm 



J 



AQ 
An ' 



(11) 



The last two co lumns of Tab [ 3] give t he range for the substructure enhancement factor B 
as obtained by IStrigari et al.l (j2007al ). The rate of dark matter pair-annihilation in the 
presence of substructures is effectively boosted by the factor B. While the enhancement 
range is uncertain both because of the theoretical modeling and because the precise value 
of the cut-off scale for the sm allest collapsed substr uctures depends critically on the particle 



da rk matter model (see e.g. iProfumo et al.l l2006l ). we refer here to the results reported 



by IStrigari et al.l (l2007al ) in their Fig. 5, and assume a small substructure cut-off, Mcut ~ 
10~^ Mq. T he results of re c ent h igh-resolution N-body simulation of Milky- Way size galaxies 
reported in iKuhlen et al.l (120081 ) indicate that the boost factor from substructures might 
be significantly smaller. For subhalo masses of the size of the dSph galaxies we consider 
here, and ext r apola ting the matter power spectrum with mass functions dN/dM oc M~", 
Kuhlen et al.l (120081 ) finds boost factors ranging from ~ 1 to ~ 40. The latter value refers 
to Mcut ~ 10^^^ Mq and to a = 2.1, and should thus be regraded as an upper bound. 
The disruption of small substructures by stellar encounters can also play a relevant role, as 
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discussed e.g. in IZhao et al.l (120071). It is beyond the scope of this analysis to assess the 



stabihty of the ranges quoted in IStrigari et ahl (l2007al ) against different N-body simulation 



results, extrapolations for the matter power spectrum at small scales and the mentioned 
particle dark matter uncertainties. We however warn the reader that the figures quoted in 
the last two columns of Tab [3] should be regraded as optimistic upper limits. 

Fixing the dark matter density profile allows us to translate the X-ray flux limits given 
in the preceding section into actual constraints on the particle dark matter models. Our 
reference choices for the diffusion setup were specified above in Sec. ^ but the crucial depen- 
dence on the diffusion parameter D will be further assessed here. Since it gives the strongest 
constraints on the X-ray flux, we choose to normalize our constraints to the Fornax dSph, 
with the reference dark matter setup specified in the second line of Tab. [21 

The left panel of Fig. [3] illustrates the relative strength of the constraints we get here 
for the three dSph under consideration, factoring in the uncertainty from the dark matter 
density distribution. A smaller number corresponds to a more stringent constraint, for 
instance on the pair-annihilation cross section for a given dark matter particle mass. The 
blue intervals correspond to the ranges on the normalization factor without substructures 
given in Tab.O The dotted red ranges indicate the expected improve ment on the constraint s 



when the effect of substructures is included, according to the model of IStrigari et al.l (l2007al ). 
The constraints obtained including substructures improve by the factors given in the last 
two columns of Tab. [31 From the figure, we deduce that the impact on dark matter models 
of X-ray observations of Fornax is comparable to that obtained from data on Ursa Minor, 
because while the latter features a larger integrated dark matter density squared by a factor 
~ 3, its flare-free XMM exposure is much lower. Data from Carina provide constraints that 
are in the same ballpark of the other two dSph, but typically less stringent by factors of a 
few. 

The right panel of Fig. [31 shows the effect on the flux of X-ray in the 0.5 to 8 keV band of 
varying the diffusion coefficient Dq in the range between 10^^ and 10^^ cm^/s. We normalize 
the flux to that obtained in our reference diffusion setup, and show the lines corresponding 
to 7 = (solid black) and 7 = 1 (red dashed). Intermediate values of 7 fall between the 
two lines. Recall that the reference diffusion setup features 7 = 0.7. The proximity of the 
two lines indicates that the precise value of 7 is much less critical to the X-ray emission 
from e^e~ produced by dark matter annihilation than the value of Do. As expected, smaller 
values of the diffusion coefficient induce a smaller loss of energetic e~^e~ , and, eventually, 
for smaller and smaller Dq the curves will converge to the value corresponding to a scenario 
where diffusion can be totally neglected. In the range we explored, diffusion can lead to a 
suppression of the signal by a factor ~ 10 for larger values of Dq, or to enhancements by 
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Table 1. Group Sample 



Dwarf 


R.A. 


Decl. 


Exposure (ks) 


Flux Limit 


Position Reference 




(J2000) 


(J2000) 


MOSl, M0S2, PN 


(photons cm~^ s~^) 




Ursa Minor 


15:09:11.3 


67:12:52 


11, 11, 7.8 


2.3E-5 


Cotton et al. 1999 


Fornax 


02:39:52.0 


-34:30:49 


67, 69, 55 


l.OE-5 


Battaglia et al. 2006 


Carina 


06:41:36.7 


-50:57:58 


19, 21, 13 


2.1E-5 


Lauberts 1982 



Note. — Flux limits are listed for the 0.5-8 keV band for an aperture of 6' radius. 



Table 2. Dark Matter Profiles 



Dwarf pf pf"" p^^^ rf rf'^ r,"^^^ D 



Ursa Minor 7.5 7.35 7.85 0.2 0.067 -0.033 66 
Fornax 7.6 7.35 7.9 0.05 0.067 -0.067 138 

Carina 7.8 7.5 8.0 -0.3 -0.23 -0.36 101 



Note. — Range and reference values for the dark mat- 
ter profiles of the three dSph galaxies under consideration. 
Columns 2, 3 and 4 indicate l ogin[/?.q/(MM kpc"^)] , columns 
5, 6 and 7 quote logj^Q[rs/kpc] (IStrigari et al.ll2007al ). and the 
last column is the reference dSph distance, in kpc. 
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Fig. 3. — Left: a comparison of the constraints on WIMP models from the three dwarf 
galaxies under consideration here, normalized to the constraint we get for Fornax with the 
central reference value for the dark matter distribution. Stronger constraints correspond to 
smaller values on the y axis. The blue solid lines indicate the ranges allowed, for an NFW 
profile, by dynamical data on the dark matter distribution of the various galaxies, assuming 
no substructures. The dotted red lines show how th ese ranges are affected accounting for dark 
matter sub-structures, according to the analysis of IStrigari et al.l (l2007al ). Given the lack of 
co nsensus on the size of the boost factor from substructres, we regard the estimate quoted 
in IStrigari et al.l (l2007al ) as an indication of the maximal enhancement from substructures 
(see the text for more details on this point). Right: the effect on the integrated X-ray flux 
over the [0.5, 8] keV range of varying the diffusion coefficient Dq, normalized to our reference 
diffusion model (green diamond). The black curve assumes 7 = 0, while the red dashed line 
7 = 1: at a given value of Dq, varying 7 between and 1 shifts the X-ray flux from the lower 
to the upper fine. 
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more than a factor 20 for smaller values of Dq compared to our reference setup. 

As alluded to above, given the dark matter density distribution, and fixing the dark mat- 
ter mass and the e~^e~ injection spectrum (via the dominant pair-annihilation final state), we 
can determine the maximal values of the pair- annihilation cross section (erf )o allowed by X- 
ray data. We carry out this exercise in Fig. H] assuming no substructure enhancement, and for 
the reference dark matter and diffusion setups, for the benchmark case of Fornax. The solid, 
dashed and dot-dashed lines correspond to the three dominant pair-annihilation final states, 
respectively bb, W^W~ and r+r~. We also show the uncertainty band (with respect to 
the bb final state) stemming from the determination of the dark matter density distribution, 
neglecting the effect of substructures, and the possibility of assuming density distribution 
profiles different from a NFW profile. Under these very conservative assumptions we are 
able to constrain interesting valu es of the pair an nihilation cross section, particularly in the 



light WIMP mass range (see e.g. iProfumd l2008bl . for a discussion of the phenomenology of 



light (rriDM <^ 1 GeV) neutralinos) . Models producing a softer e~^e~ spectrum, such as bb are 
more strongly constrained than models featuring a harder spectrum (such as r+r", where 
energetic electrons and positrons are produced in both the leptonic and in the hadronic r 
decays). Interestingly, the final states that are usually dominant in supersymmetric models, 
bb and W^W~ , give comparable constraints, well within the dark matter profile uncertain- 
ties. Since pair-annihilation modes such as r+r~ occur in supersymmetric models with a 
significant branching ratio only if, for instance, the scalar partners of the r are light, and 
are in any case subdominant with respect to the decays into heavy quarks, it is legitimate 
here to use, for the case of supersymmetry, what we obtain for a bb final state. 

We give an overview of the constraints we can place on supersymmetric models in 
Fig. [5l for different assumptions on the diffusion and the dark matter density distribution 
for the Fornax dSph galaxy. We again illustrate our constraints in the (moM, {o'v)o) param- 
eter space, and assume that the spectrum of e~^e~ is close enough to a bb final state. Dark 
shaded (orange) regions correspond to values for (mDM, {o'v)o) populated by supersymmetric 
models with a thermal relic abundance consistent with the inferred cold dark matter abun- 



dance (see e.g. iKomatsu et al.ll2008l ). Light shaded (yellow) points indicate those portions 
of the parameter space where in a standard ACDM cosmology the thermal relic abundance 
of dark matter is under-produced. In this case, we assume that either non-thermal dark 
matter production occurs, or that the Universe underwent a non-standard expansion rate 
pha se prior to Big Bang Nu cleosynthesis and in particular around the WIMP freeze-out (see 



e.g. iProfumo fc UUiol l2003l . for the discussion of one such non-standard cosmological setup 
involving a Quintessence field describing dark energy). In short: the dark shaded region 
indicates where the WIMP pair-annihilation cross section is compatible with a "vanilla" 
thermal freeze-out scenario, while the light shaded region gives models with a more extreme 
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Fig. 4. — A comparison of the constraints on the WIMP mass versus pair-annihilation 
cross section plane for different final states: bb, W^W~ and r+r^, neglecting dark matter 
substructures, and for the case of Fornax. The gray band indicates the uncertainty (for 
definiteness around the bb line) in the dark matter profile modeling. 
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Fornax 




Fig. 5. — Overview of the constraints from X-ray observations of Fornax on supersymmetric 
WIMP dark matter models. We show the hmits obtained using our conservative reference 
diffusion setup and a diffusion coefficient Dq = 10^^ cm^/s. The sohd hues correspond to the 
case of no substructures, while the dashed lines indicate the range where one could expect 
the limit to be set when substructures are included. The yellow area corresponds to values 
of the dark matter mass and annihilation cross section found for neutralino models within 
the minimal supersymmetric extension of the Standard Model. The orange area indicates 
those supersymmetric models that also produce a thermal relic neutralino abundance in 
agreement with the inferred cold dark matter density. The solid green and indigo lines 
locate the predictions for "vanilla" Wino and Higgsino dark matter models. 
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annihilation rate, compatible with the cold dark matter density provided additional assump- 
tions are made on the WIMP p r oduct ion mechanism beyond standard thermal generation. 
We refer the reader to iProfumd (120051 ) for details on the scan of the parameter space of the 
minimal supersymmetric extension of the Standard Model employed to obtain the shaded 
regions shown in Fig. [5l 

The upper solid black line reproduces the constraint shown in Fig. HJ and assumes 
no substructures and a conservative value for the diffusion coefficient (i.e. the reference 
model matching the median diffusion coefficient for the Milky Way). We compare the X-ray 
constraints with the EGRET limit on the gamma-ray flux, indicated with the nearby dot- 
dashed red line. For EGRET, we assu me a point-source sen sitivity of around 10~®cm~^s~^, 
and an angular acceptance of 1 deg (IHartman et al.lll999l ). We see that X-ray limits are 
at least as constraining as existing gamma-ray limits, even with conservative assumptions 
on the e^e~ diffusion setup. The blue solid line roughly one order of magnitude below 
the two lines discussed above indicates the constraints for a reduced diffusion coefficient, 
Dq = 10^^ cm^/s, our estimate for the diffusion coefficient scaled down from the Milky 
Way to the scale of a dSph galaxy (Sec. 2). With this diffusion setup, existing X-ray limits 
are one order of magnitude better than existing gamma-ray limits, and can be in some 
cases competitive with the GLAST anticipat ed performance (a si milar result was recently 
obtained for the case of the galactic center by lRegis &: UUiol (120081 ) . although the dominant 
X-ray production mechanism in that case is synchrotron radiation and not Inverse Compton 
scattering) . 

Adding substructures to the description of the dark matter density pushes the lim- 
its well within the area where WIMPs are produced thermally in the early universe, even 
with a conservative diffusion setup. This statement is substantiated in Fig. [5] by the two 
sets of black and blue dashed lines, corresponding, respectively for Dq = 1.1 x 10^'' cm^/s 
and Dq = 10^^ cm^/s, to the substructure boost factors indicated in Tab. O For refer- 
ence, we indicate with solid green and indigo lines the expectation for dominantly wino- 
like and higgsino-like lightest neutralinos. The case of wino-like dark matter is ubiquitous 
in so-called anomaly- mediated supersymmetry-breaking scenarios, where the SU(2) soft- 
super symmet ry breaking gaugino na ass is much smaller than its U(l) hyper-charge counter- 
part (see e.g. iMoroi &: Randallll2000l ). Wino dark matter is constrained, in the most extreme 
setup of substructure enhanced and suppressed diffusion, up to masses around 400 GeV. 
Hig gsino dark matter is also ubiquitous in several models of supersymmetry breaking (see 



Baer et al.ll2005f ). 



Jungrnan et al.lll996l ). including the focus point region of minimal supergravity (see e.g. 



An important point we wish to make here is that the constraints we show here are 
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consistent with all availab le particle physics bounds on dark matter. In particular, it was 
shown in iProfumd (l2008al ) (see the upper panel of Fig. 2) that conservative limits in the 
{rriDM, {o'v)q) plane for 10 < mDM/GeV < 1000 range between {av)o/ {cm^s~^) < 10~^^ and 
(cjf )o/(cm^s~^) < 10~^^. Given uncertainties on cosmic ray propagation in the galaxy, these 
limits were obtained assuming a cored dark matter density profile and considering gamma- 
ra y data from EGRET and H. E.S.S. on t he ga lactic center region, according to the analysis 
of ICesarini et al.l (120041 ) and of iProfumo J2005h. Li mits on the dark matter pair annihilation 
cross section were also obtained by iBeacom et al.l (120071 ) by comparing the cosmic diffuse 
neutrino signal that would result from dark matter pair annihilation and com paring it to the 



measu red terrestrial atmospheric neutrino background. The limits quoted by iBeacom et al. 
J2007h are also at the level of {av)o/{cmh-^) < 2 x lO-^i in the 10 < moM/GeV < 1000 
range. In short, the limits we quote here fall in a region of the particle dark matter parameter 
space which is not ruled out by existing data on indirect detection. 

To summarize, limits on the diffuse X-ray flux from dSph galaxies effectively constrain 
WIMP dark matter models. The strength of these constraints depends sensitively on the 
dark matter density distribution (and more specifically on the effect of dark matter sub- 
structures) and on the diffusion model. Less conservative assumptions on either one of those 
points pushes the constraints into interesting regions of the parameter space of particle dark 
matter models (see Fig. [5]). We showed that currently available X-ray and gamma-ray data 
from nearby dSph galaxies put comparable constraints on particle dark matter indirect de- 
tection even for very conservative diffusion setups; a suppression of cosmic ray diffusion 
appropriate for a dSph scale galaxy makes X-ray data more sensitive to particle dark matter 
pair annihilation by about one order of magnitude compared to gamma rays. In the next 
section we delineate a comparison between the soon to be launched gamma-ray telescope 
GLAST and future X-ray missions in terms of their sensitivity on extra-galactic particle 
dark matter searches. 



5. Summary and Discussion 

We pointed out that local dSph galaxies are an ideal environment for particle dark 
matter searches with X-rays. We described how X-rays are produced as secondary radiation 
in Inverse Compton scattering off cosmic microwave background photons of electrons and 
positrons resulting from particle dark matter annihilation. The resulting spectrum is only 
mildly dependent on the details of the particle dark matter model (the dark matter mass and 
the dominant final state into which it pair annihilates), and it is, generically, hard (spectral 
index smaller than ~ 1.5). The normalization of the emission depends on (1) the particle 
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dark matter pair annihilation rate, (2) the diffusion setup, and (3) the dark matter density 
distribution. For reasonable choices of these three a priori unknown inputs of the problem, 
the X-ray emission is potentially within reach of current X-ray detectors. Interestingly 
enough, the shape of the spectral energy distribution indicates that for dSph galaxies X-rays 
have a comparable, if not better, sensitivity to indirect dark matter detection than gamma 
rays. 

We used XMM-Newton archival data on three Local Group dSph galaxies, Ursa Minor, 
Fornax and Carina, to search for the diffuse X-ray emission expected from dark matter 
annihilation. We studied the optimal energy and radial range to search for this type of 
emission, and concluded that for XMM-Newton and for the dSph galaxies under investigation 
these correspond to an energy band between 0.5 and 8 keV and to a radius of around 6'. 
We do not find any significant signal over background, and this, in turn, was turned into 
constraints on particle dark matter models. The best constraints result from both the Fornax 
and the Ursa Minor observations, while data from Carina result in bounds that arc a factor 
of a few weaker. Ursa Minor features the largest dark matter density, making it the best 
candidate target, but has the shortest usable XMM exposure. 

In determining the impact on particle dark matter searches of our X-ray constraints, we 
pointed out the uncertainties resulting from the modeling of cosmic ray diffusion processes, 
and from the dark matter distribution. In particular, including dark matter substructures can 
boost our constraints significantly. We phrase the bounds we obtain in terms of the dominant 
dark matter annihilation final states. For those final states relevant for specific dark matter 
models, such as supersymmetry, the constraints on the mass versus pair annihilation plane are 
very similar. We then proceeded to examine how X-ray constraints on particle dark matter 
annihilation in local dSph galaxies limit the available parameter space of supersymmetric 
dark matter. In the most conservative setup only models with rather large annihilation cross 
sections are excluded. Assuming a smaller diffusion coefficient, or factoring in the effect of 
dark matter substructures, our constraints fall well within the interesting region where the 
supersymmetric dark matter can be a thermal relic from the early universe. Also, we were 
able to set limits on particular supersymmetric dark matter scenarios, such as Wino or 
Higgsino lightest neutralino dark matter. 

An important result of the present analysis is that even assuming a conservative diffusion 
setup the sensitivity of X-rays and of gamma rays to particle dark matter annihilation in 
dSph galaxies are comparable. This fact has two-fold implications: on the one hand, if longer 
observations of dSph galaxies were carried out with existing telescopes, it is possible that the 
first astronomical signature of particle dark matter annihilation would come from X-rays. 
Secondly, should a signature be detected in the future with gamma-ray telescopes, it would 
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be extremely important to confirm the nature of the signal via X-ray observations. 

In this respect, it is relevant to comment here on how future gamma-ray and X-ray 
telescopes will improve indirect dark matter searches through observations of nearby dSph 
galaxies. The LAT instrument on-board the soon to be launched GLAST satellite will extend 
the gamma-ray energy range available to EGRET, with tremendously increased effective area 
and energy as well as angular resolution. GLAST will be an ideal telescope to search for 
dark matter annihilation in dSph galaxies. Assuming a diffuse background flux of 1.5 x 10~^ 
photons cm~^ s~^ sr~^ integrated above 0.1 GeV, and an effective spectral index in the 
gamma-ray band of 2.1, we flnd that the GLAST LAT sensitivitjfl from 5 years of data will 
improve over the EGRET point-source sensitivity by large factors. In the mass versus pair 
annihilation cross section plane, and assuming a soft gamma ray spectrum (bb), GLAST will 
improve over EGRET by factors ranging between ~ 10 and ~ 100, the flrst corresponding 
to a light rriDM ~ 10 GeV dark matter particle, and the latter to a heavy one (mDM ~ 1000 
GeV). Assuming a harder gamm a-ray spectrum, as approp riate for other dark matter models 



ap 

e.g. universal extra dimensions. iHooper fc Profumdl2007l ). the GLAST performance will be 



factors between 30 and 300 better than EGRET. A signal of dark matter pair annihilations 
in gamma-rays appears therefore very promising with GLAST. If detected, such a source 
would need to be conflrmed in its nature, and using a multi-wavelength approach is one of 
the most promising strategies. 

Future X-ray telescopes, like Constellation-X and XEUS, will also have greatly increased 
effective areas with respect to current instruments. Using the currently available projections 
for the effective areas and backgrounds of these telescopeqj, we estimate that the X-ray limits 
(0.5-8 keV band) placed by a 100 ksec observation of Ursa Minor with Constellation-X or 
XEUS would improve over the limits placed in this paper by factors of roughly 35 and 70, 
respectively. Thus even for a conservative diffusion model, the future generation of X-ray 
telescopes will place similar constraints on dark matter annihilation from dwarfs as GLAST, 
stronger constraints at particle masses below a few hundred GeV. In addition, a signal from 
GLAST could be conflrmed with X-ray observations. 

Observations at other wavelengths will also be of great relevance to identify particle 
dark matter and its properties. In particular, a diffuse radio signal should also be part of the 
multiwavelength yield of particle dark matter annihilation. However, the level of the radio 
emission is crucially dependent upon the average magnetic fleld, which adds further uncer- 



^http:/ /www-glast. slac.stanford.edu/software/IS/glast Jat_perforniance.htm 
'^http://sci.esa.int/science-e/www/object/index.cfm?fobjectid=42273 



and |http : / / constellation, gsfc.nasa. gov / resources /response-matrices / index, html 



- 26 - 



tainties both in setting constraints and in understanding the nature of particle dark matter, 
should a signal be detected. Observations in the hard X-ray band would also be useful; 
however, as opposed to cluster of galaxies, where the effect of diffusion on the dark ma tter 
multi- wavelength SED is typically mild (see IColafrancesco et al.l l2006l : iProfumd l2008al ) . in 
dSph galaxies high energy electrons and positrons escape more efficiently from the diffusive 
region, suppressing the hard X-ray emission. 



Strigari et al.l (j2007bl ) recently estimated the gamma-ray flux from dark matter annihi- 



lation in newly discovered, extremely low luminosity and dark matter dominated Milky Way 
dwarf galaxies. The expected gamma-ray flux from these nearby galaxies can be larger than 
that from previously known dwarfs, depending on the so far poorly known dark matter den- 
sity distribution. Our results indicate that these objects would also potentially be excellent 
targets for particle dark matter searches with X-ray observations. 

In short, we showed that X-rays can play an important role in exploring the nature 
of particle dark matter and in pinpointing its properties. This role is complementary, but 
not subsidiary, to searches with gamma rays, and we believe very exciting results at both 
frequencies might be just around the corner. 



We acknowledge useful discussions with Fiorenza Donato and Piero Ullio on cosmic 
ray diffusion. T.E.J, is grateful for support from the Alexander F. Morrison Fellowship, 
administered through the University of California Observatories and the Regents of the 
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Table 3. Dark Matter Density Squared Integrals and Substructure Enhancement Ranges 



Dwarf Jrci Jmin <-^max -^Iqw -^high 



Ursa Minor 3.1 1.1 6.7 25 89 
Fornax 1.2 0.42 2.49 50 159 

Carina 0.53 0.24 1.13 50 80 



Note. — The hne-of-sight integral of dark 
matter density squared, see Eq. ffTOl) . in units 
of 10^^ GeV^ cm~^, and the range for the 
substructure boost factor B as estimated in 



Strigari et al.l (l2007al ) , for the three dSph galax- 



ies under consideration here. 
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